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ABSTRACT 


^~rhe  response  of  the  seasonal  thermocline  formation  and 
miring  to  prescribed  vertical  interior  motion  is  examined. 
For  the  annual  and  shorter  period  interior  motion  cases 
studied,  the  response  was  strongest  for  the  annual  period. 
For  an  oscillatory  vertical  motion  having  a  15  m  amplitude 
at  a  75  m  depth,  warm  sea  surface  temperature  anomalies  of 
up  to  2.  16?C  and  cold  anomalies  of  up  to  2.  04^  are  found, 
depending  on  the  phasa  difference  between  the  interior 
motion  and  the  annual  heat  cycle.  The  net  phase-averaged 
effect  on  the  mixed  layer  for  annual  period  interior  motion 
is  a  reduction  in  vertical  mixing.  Higher  frequency  motion 
produces  a  net  enhancement  of  the  mixing. 
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I.  INTRODUCTION 


Interior  vertical  motion  can  play  a  significant  role  in 
the  mechanical  energy  budget  of  the  ocean  boundary  layer. 
Linden  (1975)  carried  out  experiments  to  observe  mixed  layer 
deepening  in  a  stratified  fluid.  He  found  that  internal 
waves  are  important  in  the  energetics  of  the  mixed  layer  and 
that  the  mixing  process  may  provide  a  significant  source  of 
energy  for  the  generation  of  internal  waves.  In  this  thesis 
the  interior  motion  is  prescribed  and  not  affected  by  the 
mixed  layer  processes.  However,  the  modelling  of  the  inter¬ 
action  of  free  internal  waves  with  mixed  layer  processes  is 
an  important  area  for  future  study.  Linden  points  cut  that 
he  did  not  account  for  waves  generated  elsewhere  which 
propagate  into  the  region  of  consideration.  Also,  the 
internal  waves  considered  by  Linden  were  cf  relatively  high 
frequency. 

De  Szceke  (1980)  found  wind  stress  curl- driven  vertical 
advecticn  to  be  as  important  as  the  entrainment  rate  in 
mixed  layer  dynamics.  Cus hman-Hcisin  (1981)  likewise  found 
wind-driven  advecticn  to  be  as  important  as  wind-driver, 
mixing  in  mixed  layer  dynamics  and  temperature  front 
formation. 

Much  work  has  been  done  to  study  the  ocean  response  to 
storms.  Price  et  al.  (1978),  in  their  work  or  the  mixed 
layer  response  to  stcrms,  found  it  necessary  to  include 
vertical  advection.  Price  (1981)  found  that  upwellicg  will 
significantly  enhance  the  sea-surface  temperature  response 
to  a  slowly  moving  hurricane.  Adamec  at  al.  (1981),  using  a 
two-dimensional  model  to  study  the  oceanic  response  to  a 
hurricane,  found  a  strong  interdependence  between  mixing  and 
advection.  Sreatbatch  (1  983)  also  found  advecticn  to  be 


important  in  increasing  maximum  cooling  due  no  the  passage 
of  a  storm. 

Only  recently  has  much  been  done  to  include  long  period 
interior  motion  in  mixed  layer  models.  Stevenson 
(1980,1981)  examined  the  response  of  the  mixed  layer  when 
perturbed  by  linear  Bossby  waves.  He  found  no  mixed  layer 
response  during  the  first  half  of  the  heating  season  while 
the  mixed  layer  is  shallowing  and  a  maximum  sea  surface 
temperature  response  at  the  end  of  the  heating  season  and 
beginning  of  the  cooling  season.  Stevenson  found  that  the 
maximum  response  of  mixed  layer  depth  to  the  wave  occurred 
at  or  near  the  end  of  the  cooling  season.  Also,  he  found 
that  the  lowest  frequency  waves  produced  the  greatest  mixed 
layer  depth  and  temperature  responses.  This  study  differs 
from  that  of  Stevenson  in  that  a  finite  amplitude  prescribed 
interior  motion  is  used  rather  than  a  linear  Rossby  wave. 
This  produces  a  phase  dependent  asymmetry  in  the  mixed  layer 
response  to  the  interior  motion. 

Burger  (1982)  successfully  introduced  tidal  period 
vertical  advection  into  the  Garwood  (1976,1977)  one¬ 
dimensional  mixed  layer  model.  He  found  that  the  interior 
motion  had  a  significant  effect  on  the  near-surface  tempera¬ 
ture  field  and  that  the  inclusion  of  the  interior  motion  can 
improve  the  model  for  single  station  forecasting. 

In  this  study,  the  Garwood  (1976,1977)  one-dimensional 
ocean  planetary  boundary  layer  model  is  used.  Vertical 

advection  was  added  to  the  model.  The  annual  heat  cycle  is 
a  sine  function  and  wind  stress  is  held  constant  at  one  dyne 
per  square  centimeter.  By  keeping  the  surface  boundary 

conditions  simple,  the  effect  of  the  interior  motior  can  be 
more  readily  observed. 

The  interior  vertical  motion  is  a  prescribed,  forced 
interior  wave  which  is  linear  with  depth  and  sinusoidal  in 
time  : 
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(1.1) 


w  (1  ,  t  ) 


~  w,  s.h 

D 


mi rc 


where  D  is  a  reference  depth,  200  meters,  and  w0  is  the 
amplitude  of  the  prescribed  vertical  motion  at  depth,  D. 
The  periods,  T,  studied  are  one  year,  one-half  year,  one- 
third  year  and  one-quarter  year.  The  phase  difference,  <t>  , 
between  the  interior  motion  and  the  annual  heating  cycle  is 
varied  from  0°  to  360°  at  increments  of  22.5°.  This 
interior  motion  is  general  and  may  be  representative  of  a 
number  cf  possible  long  period  phenomena,  for  example, 
internal  Kcssby  waves  or  perhaps  Elcman  pumping  on  seasonal 
scales. 

The  oceanic  variables  of  interest  are  mixed  layer  depth, 
h,  sea  surface  temperature  which  equals  mixed  layer  tempera¬ 
ture,  fMl_,  the  heat  content,  H,  and  the  potential  energy,  P. 
These  variables  are  contoured  as  functions  of  time  and  phase 
for  each  cf  the  pericds  studied.  Also,  contours  are  made  of 
the  differences  of  these  variables  relative  to  the  case  of 
no  interior  motion.  These  contours  are  used  to  observe  the 
effect  of  the  interior  motion  on  the  mixed  layer  processes. 

Furthermore,  the  phase,  <p ,  in  equation  (1.1)  can  be 
replaced  by  the  expression: 


3  TT  x 

L. 


(1.2) 


uh9re  x  is  horizontal  displacement  and  L  is  some  typical 
wavelength.  Then,  at  any  fixed  x,  equation  (1.2)  represents 
the  vertical  motion  due  to  the  passage  cf  a  travelling  wave. 
The  wavelength,  L,  might  be,  for  example,  300  km  for  an 
annual  baroclinic  Bossby  wave  in  the  central  north  Pacific 
(Kang  and  Hagaard,  1980).  The  phase  in  equation  (1.1)  can 
thus  be  related  to  horizontal  location.  The  contours  of  h, 
fML,  H  and  P  versus  phase  and  time  then  indicate  horizontal 
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structure  in  these  variables  due  to  rhe  presence  of 
wave. 
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II.  jJODEi  AND  BgOATIDIS 


A.  TBS  MODEL 

1 .  Sobs  general  Considerations 

For  this  study  the  Garwood  (1976,  1977)  one- 
dimensional  bulk  mixed  layer  modal  is  used.  The  model  is 
based  on  the  Navier-Stokes  equations  of  motion,  assumes  the 
continuity  equation  for  an  incompressible  fluid  and  makes 
the  one-dimensicnal  approximation.  The  Garwood  model 

differs  from  previous  bulk  mixed  layer  models  in  several 
features.  The  time  scale  for  viscous  dissipation  cf  turbu¬ 
lent  kinetic  energy,  THE,  is  dependant  on  planetary  rota¬ 
tion.  This  allows  for  a  cyclical  steady  state  on  an  annual 
basis  by  increasing  dissipation  for  deeper  mixed  layers. 
Both  downward  and  upward  turbulent  heat  fluxes  are  important 
ir.  this  model.  Furthermore,  entrainment  depends  on  the 
relative  distribution  of  TKE  between  vertical  and  horizontal 
components. 

Before  continuing,  a  comment  on  notation  is  in 
order.  An  overbar  cn  a  variable  will  indicate  the  horizon¬ 
tally  averaged  mean.  A  prime  on  a  variable  will  indicate 
the  fluctuation  about  this  horizontal  mean.  Angled 
brackets,  <  >,  will  indicate  vertically  averaged,  mixed 

layer  mean. 

2.  Equation  of  State 

In  the  Garwood  modal,  the  equation  of  state  has  the 

form : 
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p=  e*  [ '  -  -(t-t.)  +  p(s-  s.)] 


(2.1) 


The  thermal  expansion  coefficient,  ^  f  is  a  function  of  T 
and  is  re-computed  daily.  The  constants  p,,  T.  and  S0  are 
typical  density,  temperature  and  salinity  values,  respec¬ 
tively.  In  this  study,  salinity  is  fixed  (precipitation  and 
evaporation  are  disallowed)  and  set  equal  to  the  value  of 
S, .  Equation  (2.1)  then  becomes 

p  =  /»*[*  -  * C1" ' T° .  (2-2) 


B.  THE  HEAT  EQUATION 

i  •  i filial 

The  time  rate  cf  change  of  temperature  is  given  by 
the  one-dimensional  heat  equation: 


2T  _ 

-  -1 

i-  | 

[4-Vwix.K  . 

(2.3) 

Jt  ~ 

^ ecr } 

The  first  term  on  the  right  is  vertical  turbulent  heat  flux, 
the  second  is  solar  insolation,  the  third  is  vertical 
temperature  advection  and  the  last,  diffusion.  All  of  the 
solar  insolation  will  be  assumed  to  be  absorbed  at  the 
surface  (in  the  first  few  centimeters)  of  the  ocean.  Hence, 
the  solar  insolation  term  will  be  dropped  in  favor  of  a 
surface  boundary  condition  on  the  turbulent  heat  flux. 
Equation  (2.3)  then  becomes 


(2.4) 
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Ir.  this  model  ,  mixed  layer  temperature  is  assumed 
to  be  vertically  homogeneous.  Equation  (2.4)  then  reduces 
to 


(2.5) 


for  mixed  layer  temperature,  TMt.  Obviously,  TH4_  is  equal 
to  sea  surface  temperature  because  of  the  vertical  homoge¬ 
neity  assumption  in  the  mixed  layer. 

The  surface  boundary  condition  on  heat  flux  is  given 
by: 


(2.6) 


where  the  zero  subscript  indicates  the  value  of  w  T'  at  z=0. 
The  net  surface  heat  flux,  Q/j°.or,  is  prescribed  to  be 
harmonic  in  time: 


(2.7) 


The  amplitude  of  this  sine  function  is  set  equal  to  a 
representative  mid-latitude  value  of  h'lQ-s'c  cm/sec.  In 
this  study,  heat  flux  is  positive  downward.  Time  t=0  is  the 
start  of  the  '’heating  season",  when  Q  first  becomes 
positive . 

The  boundary  condition  on  heat  flux  at  the  base  of 
the  mixed  layer  is  given  by  the  entrainment  heat  flux. 
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(2.8 


tw'T*  )K  =  -  AT  w*  , 

where  the  subscript,  h,  indicates  the  value  of  w'  T'  at  z=-h , 
the  base  of  the  nixed  layer.  The  variable.  We  is  the 
entrainment  rate  and  dT  is  the  temperature  jump  across  the 
entrainment  zone  at  the  base  of  the  mixed  layer.  The  turbu¬ 
lent  heat  flux  is  linear  with  depth  in  the  mixed  layer  (due 
to  the  assumption  of  a  well-mixed  layer),  mcnct cnically 
decreases  with  depth  in  the  entrainment  zone  and  is  zero 
below.  The  profile  of  turbulent  heat  flux  is  shewn  in 
Figure  2.1.  The  net  surface  heat  flux  varies  from  positive 


Figure  2.1  Turbulent  Heat  Flux  Profile. 

to  negative  values  according  to  eguation  (2.7) .  The 
entrainment  heat  flux  term,  -£T  We,  is  always  negative  or 
zero.  This  i3  because  the  entrainment  rats  is  positive  when 
the  mixed  layer  is  deepening  but  "turns  off"  and  is  zero  for 
non-entraining  shallowing  cases.  The  thickness  of  the 


8 


entrainaert  zone,  6,  is  assumed  to  be  much  less  than  the 
depth  of  the  aixed  layer,  h. 

5)  now  becomes 

Q  ATWc  (2.9) 

/°«  CP  h  K 

Note  that  the  advecticn  tar  a  is  absent.  Thus,  the  interior 
aotion  does  not  directly  affect  the  aixed  layer  temperature. 
Instead,  tha  interior  aotion  affects  the  aixed  layer 
temperature  indirectly  through  the  variables  h,  at  and  We. 


Equation  (2. 


dT, 


<i  u 


Z  t 


Figure  2.2  Initial  Temperature  Profile 


The  initial  temperature  profile  is  given  in  Figure 
2.2.  The  reference  depth,  D,  is  203  at  and  will  be  used  in 
sobs  of  the  equations  in  later  sections. 

2.  Scale  Analysis  an£  Siapllf  ication  of  the  Heat 

Hguation 

Consider  now  the  relative  magnitudes  of  the  terms  of 
equation  (2.4).  In  the  mixed  layer,  the  temperature  may 
change  by  as  much  as  10  °C  in  a  half  year.  The  time  rate  of 
change  of  mixed  layer  temperature  would  then  be  or.  the  order 
of  10"*  °C/sec  or  less.  3elow  the  mixed  layer  the  time  rate 
of  change  of  temperature  of  fluid  particles  will  generally 
be  smaller. 

The  first  term  on  the  right  of  equation  (2.4)  is  the 
turbulent  heat  flux  term.  In  the  mixed  layer,  the  only  term 
on  the  right  side  of  equation  (2.5)  is  the  turbulent  heat 
flux.  Thus  this  term  will  scale  with  the  seasonal  time  rate 
of  change  of  mixed  layer  temperature,  10'*  Cc/sec.  The 
prescribed  surface  heat  flux  boundary  condition,  Q/ p£ cp  ,  has 
extreme  values  of  t4*10'3<,C  cm/sec.  The  value  of  h  varies 
from  around  30  a  to  around  100  m.  Thus,  the  first  term  on 
the  right  side  cf  equation  (2.9)  has  a  range  of  values 
between  ±10'6*C/sec.  The  second  term  on  the  right  of  equa¬ 
tion  (2.9)  has  roughly  the  same  order  of  magnitude,  but  it 
can  only  be  negative.  Below  the  mixed  layer  and  entrainment 
tone,  that  is,  below  the  flux  region,  the  turbulent  heat 
flux  terir  vanishes. 

next  consider  the  advection  term  in  equation  (2.4). 
Recall  that  this  term  is  zero  in  the  mixed  layer.  The 
gradient  of  the  temperature  profile  just  below  the  mixed 
layer  is  initially  about  10'*°C/cm.  During  the  heating 
season  the  temperature  at  the  top  of  the  thermocline  will 
increase,  but  the  vertical  extent  of  the  thermocline  will 
also  increase  as  the  mixed  layer  shallows.  This  term  should 
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thus  te  expected  to  have  the  sane  order  of  magnitude.  The 
amplitude  of  W  will  be  prescribed  according  to  equation 
(2.19)  in  section  D.  The  maximum  value  at  200  m  depends  on 
the  wave  or  oscillation  frequency,  but  it  will  be  on  the 
order  of  10-3  cm/sec.  This  amplitude  will  be  smaller  at 
lesser  depths.  The  result  is  that  the  advection  term  of 
equation  (2.4)  should  be  of  the  order  of  10"‘°C/sec  below 
the  mixed  layer. 

Finally,  consider  the  diffusion  term  of  equation 
(2.4).  This  term  vanishes  in  the  mixed  layer.  Below  the 
mixed  layer  a  reasonable  scale  for  the  second  derivative  of 
temperature  is  on  the  order  of  4cC/(5000  cm)2.  A  typical 
range  of  values  for  the  eddy  diffusivity,  «  ,  is  0.1  to 
0.5  cm2/sec.  Thus,  the  diffusion  term  has  a  magnitude  of 
about  10-8  to  10_7°C/sec.  The  diffusion  term  is  much 
smaller  than  the  other  terms  of  equation  (2.4)  and  will, 
from  now  on,  be  neglected.  The  heat  equation  (2.4)  thus  can 
te  simplified  to 


C.  MIXED  LAYER  DEPTH 

In  the  model  of  Garwood  (1977),  the  mixed  layer  either 
shallows  to  a  retreat  depth,  hr,  or  deepens  at  the  entrain¬ 
ment  rate.  We.  The  retreat  dapth  is  determined  by  assuming 
steady  state  conditions  for  the  turbulent  kinetic  energy 
budget,  letting  entrainment  heat  flux  vanish  and  letting 
<i Tr2>  approach  zero.  The  resulting  equation  for  retreat 
depth  is 
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(2.11) 


(Garvcod,  1977).  In  this  aquation,  m,  and  i,  are  dimension- 
lass  model  constants,  <x  is  the  thermal  expansion  coefficient 
and  <E>  is  the  mixed  layer  mean  turbulent  kinetic  energy. 
Also,  u„,  -.he  friction  velocity,  is  the  square  root  of  the 
surface  wind  stress  divided  by  density.  Throughout  this 
study,  it  is  set  equal  to  one  centimeter  per  second.  Since 
mixed  layer  mean  turbulent  kinetic  energy  depends  or.  surface 
boundary  conditions,  the  retreat  depth  is  independent  of  the 
prescribed  interior  action. 

In  the  case  of  active  entrainment,  the  mixed  layer  deep¬ 
ening  rate  in  the  model  depends  on  the  entrainment  rate  and 
WK,  the  mean  vertical  velocity  at  the  bottom  of  the  mixed 
layer  due  to  the  prescribed  interior  motion: 


ih_ 
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(2.  12) 


The  entrainment  rate.  We,  depends  on  the  entrainment  zone 
temperature  jump,  £T,  and  the  entrainment  heat  flux  (see 
equation  (2.8)).  The  entrainment  heat  flux  is  modelled  by 
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(Garwood,  1977),  where  m4  is  a  model 
entrainment  rate  can  then  be  expressed  as 
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Therefore,  equations  (2.12)  and  (2.14)  give  a  solution  for 
Jh/<)t  if  Wh,  </  z>  and  <E>  are  known. 

For  moderately  strong  wind  stress,  mixed  layer  mean 
turbulent  kinetic  energy,  <1>,  depends  primarily  on  surface 
wind  stress,  which  is  held  constant  here.  The  vertical 
component  of  turbulent  kinetic  energy,  lyr~2)l/'2,  depends 
primarily  cn  surface  heat  flux  and  surface  wind  stress. 
When  the  surface  heat  flux  becomes  large  enough,  this  term 
becomes  zero,  stopping  tha  entrainment  process  and  permit¬ 
ting  shallowing.  These  two  terms  are  unaffected  by  the 
interior  motion.  The  entrainment  rate  is  affected  by  the 
mean  vertical  motion  indirectly  as  a  result  of  changes  in 
mixed  layer  depth,  h,  and  entrainment  zone  temperature  jump, 
AT. 


D.  IBTEBIOB  VERTICAL  HOTION 


Consider  for  the  moment  a  two  layer,  deep  ocean  with  the 
upper  layer  about  200  m  deep  and  the  lower  layer  about 
5000  a  deep.  Assuming  the  fluid  is  incompressible,  the 
continuity  equation  is  differentiated  with  respect  to  z: 


J  *  *  4  d  y  <#  •£ 


(2.  15) 


where  0,  V  and  W  are  the  mean  velocity  components  in  either 
layer.  Furthermore,  assume  that  horizontal  momentum  is 
vertically  well  mixed  in  the  upper  layer.  Thus, 


(2.  16) 
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Than  V  must  ba  linear  with  depth  in  tne  upper  layer.  With  a 
rigid-lid  approximation,  W  *  0  at  the  surface. 

Our  interest  is  confined  to  the  upper  200  m  and  a  linear 
vertical  profile  for  w  is  a  good  approximation  for  long 
period  interior  motion.  The  functional  dependence  of  w  on  z 


(z)  -  - 


(2.17) 


where  D  =  200  m  is  the  reference  depth  and  W^  is  the 
interior  vertical  motion  at  200  a  given  by  a  sinusiodal 
function  of  time,  namely 
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The  amplitude,  wa  ,  is  defined  by 
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motion. 


For  the  one-year  period  interior  motion, 

W0  =  8  *  10* 4  cm/sac  *  69  cm/day.  The  value  of  Wc  will  be  2, 
3  and  4  times  this  for  the  one-half,  one-third  and  one- 
quarter  year  periods,  respectively.  These  values  of  w0  are 
selected  so  that  the  displacement  amplitude  at  75  m  is 
approximately  15  m.  The  actual  values  of  the  displacement 
amplitude  depend  on  phase  and  will  be  discussed  below. 

Combining  equations  (2.17)  and  (2.18)  yields  the 
equation  fcr  the  prescribed  interior  motion. 


The  formula  for  the  particle  position  due  to  the  interior 
motion  can  be  obtained  by  noting  that  W  -  dz/dt,  ar.i  is 
given  by 
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For  any  fixed  phase  the  maximum  and  minimum  values  of  z  are 
given  by 
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The  maximum  peak-to-peak  displacement  amplitude  will  occur 
when  cos  <f>  =+1  and  the  minimum  when  cos  <p  =- 1 .  For  the 
vertical  motion  amplitude  given  by  equation  (2.19),  a 
particle  with  an  initial  position  at  z=-75  m  will  experience 
a  displacement  amplitude  (half  of  the  peak-to-peak  value) 
between  12.4  and  18.4  m. 

The  model  time  step  is  one  hour.  This  is  much  smaller 
than  the  cne-year  or  even  one-quarter  year  periods  consid¬ 
ered.  Hence,  since  HD  does  not  vary  significantly  over  a 
one  hour  time  step,  the  integral  of  (2.  17)  can  be  approxi¬ 
mated  by 
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where  is  the  tine  increment  (cne  hour)  and  Wn  is 

obtained  from  equation  (2.18).  In  the  model,  vertical 
advection  is  handled  by  compressing  or  expanding  the 
vertical  grid  according  to  equation  (2.25).  No  significant 
difference  in  the  results  was  found  between  model  runs  using 
equation  (2.25)  and  equation  (2.21). 

The  vertical  increment  is  initially  one  meter,  but  this 
varies  as  a  result  of  the  compression  or  expansion  of  the 
vertical  grid.  This  variation  in  the  vertical  increment  is 
never  more  than  about  25  3(.  A  consequence  of  this  scheme  is 
that  the  lowest  grid  point  is  not  always  at  a  depth  of  199  m 
as  in  the  initial  profile,  Figure  2.2.  Care  must  be  taken, 
therefore,  when  performing  vertical  integration  from  the 
surface  to  the  reference  depth,  D  =  200  ra.  The  integrals 
that  will  be  calculated  in  the  next  section  involve  the 
temperature  difference,  f  -T0  .  Problems  with  varying  the 
grid  size  can  easily  be  avoided  by  setting  the  reference 
temperature,  Tc ,  equal  to  T0,  the  initial  temperature  at 
200  m.  This  way,  T-TD  will  always  be  zero  at  the  base  of 
the  vertical  grid.  The  grid  is  initially  made  deep  enough 
so  that  the  seasonal  thermocline  never  reaches  the  bottom  of 
the  grid;  hence,  T D  always  equals  its  initial  value,  6.08JC. 

Subtracting  the  reference  temperature,  f D ,  from  tempera¬ 
ture,  ?,  in  the  integrals  in  the  next  section  amounts  to 
subtracting  the  terms; 

r°  -  = 
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Since  TD  always  remains  constant,  both  expressions,  (2.26) 
and  (2.27),  will  always  remain  constant.  Subtracting 
from  T  will  amount  to  subtracting  constants  from  the 
integrals  ir.  the  next  section. 


E.  POTENTIAL  EHEBGT  AND  NEAT 

The  potential  energy  per  unit  area  for  a  column  of  water 
from  the  surface  to  depth,  D,  is  given  by 
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Employing  equation  (2.2)  with  T0  as  the  reference 
temperature,  allows  aquation  (2.28)  t  o  be  written  as 
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The  constant  in  equation  (2.29)  is 
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The  integral  in  equation  (2.29)  will  be  used  for  the  compar¬ 
ison  of  the  potential  energy  with  and  without  the  interior 
motion.  It  will  be  represented  by  the  variable  P  with  units 
of  *C  a2  where 

P  =  -  f  (  T  -  T, l)2;li  .  (2.31) 
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Similarly,  heat  content  will  be  indicated  by  the  variable  H 
with  units  °C  m  where 


H  =  j  (T-T,) 

-D 


(2.32) 


The  analysis  of  the  physical  effects  of  the  interior 
motion  can  be  facilitated  by  developing  the  expressions  for 
the  time  rate  of  change  of  H  and  of  P.  First, 
differentiating  equation  (2.32)  with  respect  to  time  gives 


o>  H 
o  c 


(2.33) 


Equation  (2.33)  can  be  written  as 


(2.34) 


Substituting  equation  (2.10)  for  th-  time  rate  of  change  of 
temperature  into  equation  (2.34)  yields 


(2.  35) 


The  first  term  on  the  right,  when  integrated,  is  just  the 
net  surface  heat  flux.  Substituting  equation  (2.17)  into 
the  second  term  on  the  right  of  (2.35)  and  integrating  by 
parts  gives 
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Using  equation  (2.26)  this  can  be  rewritten  as 
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Equation  (2.35)  can  new  be  written  as 
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The  formulation  for  the  time  rate  of  change  of  the 
potential  energy  is  similar  and  is  gi van  by 
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Substituting  equation  (2.10)  into  (2.39)  yields 


o'  *-  t\  '  ~  rs 


...  *T 


1  W  — —  c 

■  V  i 


(2.  40) 


The  first  term  in  equation  (2.4Q>  can  be  integrated  by 
parts : 


j  ?  cU  =  2(w-t'  )\  -  J  (C3"^F)  cU  .  (2.41) 
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The  first  term  on  the  right  of  aquation  (2.41)  vanishes 
because  w'  T '  *  0  at  z  3  -D.  Since  V  T'  is  zero  below  the 
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entrainment  zone,  equation  (2.41)  reduces  to 
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This  integral  may  be  broken  up  over  the  mixed  layer  and  the 
entrain  men*’  zone,  producing 
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where  < w 7 T'  >  is  the  vertically  averaged 
the  turbulent  heat  flux.  Since  wT'  is 
to  z  in  the  mixed  layer,  this  mixed 
written 


mixed  layer  mean  of 
linear  with  respect 
layer  mean  can  be 
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Consider  the  second  term  on  the  right  side  of  equation 
(2.40).  Substituting  equation  (2.17)  for  1  and  integrating 
by  parts  yields 
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using  equation  (2.27)  this  can  be  rewritten  as 

j  *  w  \  P.  <2*w6> 
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On  combining  equations  (2.40),  (2.43),  (2.44)  and  (2.46), 

the  time  rate  of  change  of  the  potential  energy  variable  now 
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becomes 
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Because  the  turbulent  heat  flux  is  monotonic  across  the 
entrainment  zone,  as  illustrated  in  Figure  2.  1: 
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This  indicates  that  the  integral  term  in  equation  (2.47) 
should  generally  be  smaller  in  magnitude  than  the  first  term 
cn  the  right  side,  unless  of  course,  the  surface  and 
entrainment  heat  fluxes  nearly  cancel. 

Actual  values  of  the  terms  of  equation  (2.47)  were 
calculated  for  one-year  period  interior  motion  at  two 
different  phases  and  the  second  term  on  the  right  of  the 
equation  was  smaller  than  the  first,  usually  by  one  order  of 
■aqnituie.  Furthermore,  the  third  term  on  the  right  was 
generally  an  order  of  magnitude  larger  than  either  of  the 
other  two.  Equation  (2.47)  will  thus  be  approximated  by 
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III.  EXPERIMENTS  AND  fiESOLTS 

A.  SUMMARY  OF  PERTINENT  EQUATIONS 

The  equations  for  time  rate  of  change  of  mixed 
depth,  mixed  layer  (and  sea  surface)  reaperature ,  hea 
potential  energy  are 
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where  (w'T')^  has  been  replaced  by  equation  (2.6) 
The  entrainment  rate  is  given  Dy 


AT  W, 
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The  prescribed  net  surface  heat  flux  is  given  by 
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and  is  plotted  as  a  function  of  time  in  Figure  3.1. 


vi/  t,’>  ^0<“ 


Hfc  TIN  G 


t  (  ii  wv  V  > 


Figure  3.1  Net  Surface  Heat  Flux. 

The  prescribed  vertical  motion  at  the  reference  depth, 
D=2 0 0  a,  is 


0  . 


(3.8) 


where  H0  is  defined  by  equation  (2.  19)  .  By  substituting 
z*-h  into  equation  (2.17),  it  can  be  seen  that  WK  varies 
directly  as  h  and  Hs. 

Some  additional  notation  will  be  useful  for  the  discus¬ 
sion  of  the  results.  To  observe  the  effect  of  the  interior 


notion,  the  values  of  h,  rMt,  h  and  ?  are  coapared  to  these 
without  interior  action.  Contours  are  plotted  of  the 
difference  of  nixed  layer  depth  with  and  without  the 
interior  notion,  likewise  for  aixed  layer  tenperature,  heat 
and  potential  energy.  Therefore,  the  following  new 
variables  are  defined: 


A  k 


-  k  - 


(3.9) 


ATml  = 


r  -  t 

hl 


Ml  -  .  , 
^  1  W  C  i 


(3.  10) 


A  H  -  H  -  W 


1C 


(3.  11) 


Ap  -  9  -  F  _ 


(3.  12) 


B.  BES01TS  WITH  MO  INTERIOR  MOTION 

The  aixed  layer  depth  has  bean  plotted  in  Figure  A. 4. 
Note  that  the  vertical  axis  scale  is  inverted  for  clarity. 
The  year  was  started  at  the  beginning  of  the  heating  season 
and  the  aixed  layer  retreats  for  the  first  quarter  year. 
Then  the  aixed  layer  deepens  slowly  for  the  next  half  year 
and  more  rapidly  in  the  last  quarter.  In  the  last  few  days 
the  nixed  layer  shallows  to  the  t=0,1  y ear , . . .  ,a-years 
cyclical  steady  state  value  of  about  75  n. 


The  mixed  layer  temperature,  f^L,  is  plotted  ir.  Figure 
A. 5.  As  expected,  the  mixed  layer  temperature  increases 
during  the  heating  season  and  decreases  during  the  coding 
season,  for  the  most  part.  However,  note  that  the  peak 
mixed  layer  temperature  occurs  around  day  170  and  T„u  smarts 
to  decrease  ten  days  before  the  end  of  the  heating  season. 
This  is  because  the  entrainment  heat  flux  term  has  begun  to 
dominate  eguation  (3.2).  Seme  entrainment  had  already 
started  fifty  days  earlier.  3y  day  170,  towards  the  end  of 
the  heating  season,  AT  has  reaohed  its  largest  value, 
7.8  “C.  Cn  day  170  the  cooling  of  the  mixed  layer  due  to 
entrainment  has  reached  0.  021  aC  per  day  while  the  surface 
heating  has  droppad  to  0.030  3C  per  day. 

The  heat  content,  plotted  in  Figure  A. 6,  depends  only  on 
surface  heat  flux  for  the  case  of  no  interior  motion  (see 
equation  (3.3)).  Any  heat  lost  from  the  mixed  layer  due  to 
entrainment  will  be  gained  by  the  region  below.  Since  heat 
is  calculated  from  the  surface  to  200  m,  entrainment  will 
have  no  effect  on  its  value. 

Potential  energy  per  unit  area  is  plotted  in  Figure  A. 7. 
This,  toe,  generally  increases  with  heating  and  decreases 
with  cooling.  However,  note  that  the  peak  value  occurs  at 
about  day  305,  25  days  into  the  cooling  season.  This  is 
because  potential  energy  is  also  increased  by  mixing  as 
indicated  by  equation  (3.4)  .  The  entrainment  term  dominates 
equation  (3.4)  from  day  170  until  day  305,  when  the  net 
surface  heat  flux  is  near  zero.  A  close  inspection  of 
Figure  A. 7  reveal3  that  the  potential  energy  at  the  end  of 
the  year  is  greater  than  the  initial  value.  Specifically, 
the  initial  vaiue  of  E  is  19608  ;C  m*;  the  value  at  the  end 
of  the  first  modal  year  is  20204 °c  m*  and  at  the  end  of  the 
second  year  is  20579  *C  m*.  This  increase  in  the  potential 
enargy  is  due  to  the  adjustment  cf  the  temperature  profile 
as  the  ocean  approaches  cyclical  steady  state. 
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Figure  A .  1  is  a  plot  cf  temperature  versus  depth  and 
time.  The  top  of  the  graph  where  the  temperature  contours 
are  vertical  is  the  nixed  layer.  The  mixed  layer  structure 
is  quite  evident  after  the  first  quarter  year  due  to  the 
strong  temperature  gradient  at  the  base  of  the  mixed  layer. 
Note  that  the  temperature  profile  at  any  depth  below  the 
mixed  layer  is  left  unchanged  until  mixing  reaches  that 
depth . 

C.  08E-YEAB  PERIOD  INTERIOR  NOTION 

The  zero  contours  of  the  prescribed  interior  motion  have 
teen  drawn  in  Figure  A. 8  with  upwelling  and  downwelling 
indicated.  The  dotted  contours  indicate  when  the  magnitude 
of  the  interior  motion  drops  to  half  of  the  maximum  value. 

During  ~he  first  quarter  year  the  mixed  layer  is  shal¬ 
lowing.  Since  the  retreat  depth  depends  primarily  on 
surface  tcundary  conditions,  mixed  layer  depth  will  be  unaf¬ 
fected  by  interior  motion  at  this  time.  This  is  borne  out 
by  the  contours  of  ih  (defined  in  equation  (3.9)),  Figure 
A.  10.  The  first  quarter  year  shows  no  difference  in  mixed 
layer  depth.  With  little  or  no  effect  on  mixed  layer  depth 
during  the  first  quarter  year,  there  will  be  little  cr  no 
difference  in  the  value  of  the  first  term  of  equation  (3.2) 
with  and  without  the  interior  motion.  The  second  term  of 
equation  (3.2)  will  be  zero  at  this  time  since  entrainment 
stops  during  shallowing.  Hence,  mixed  layer  temperature 
will  alsc  remain  unaffected  by  the  interior  motion  during 
the  first  quarter  year.  This  can  be  seen  in  Figure  A. 12. 
Heat  and  potential  energy,  on  the  other  hand,  are  influenced 
by  the  advection  of  temperature  in  the  thermocline  as  indi- 

V 

cated  by  the  presence  of  the  last  term  in  equation  (3.3)  and 
equation  (3.4)  .  Thus,  the  interior  motion  will  influence 
heat  content  and  potential  energy  even  during  the  first 
quarter  year  (Figures  A.  16  and  A.  19). 


*r  ■■  ■  —  ■■  - -  - - - 
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1  •  Effect  on  M  ixed  La^er  De p-.h 

Equation  (3.1)  indicates  rhat  downwelling  will 
enhance,  deepening,  and  upwelling  will  reduce  deepening.  In 
Figure  A. 9,  it  can  be  seen  -hat  deepening  still  dees  occur 
after  the  first  quarter  y*ar  as  in  did  without  any  interior 
motion.  However,  the  rate  of  deepening  depends  on  phase. 

At  the  90"  phase,  there  is  upwellir.g  during  the 
middle  half  of  the  model  year  (see  Figure  A.  2)  reducing  the 
rate  of  deepening.  At  day  274  down  welling  begins  and  the 
mixed  layer  begins  to  deepen  more  rapidly.  The  very  rapid 
deepening  in  the  last  quarter  of  tne  year  for  this  phase 
results  from  both  entrainment  and  downwelling  working 
together.  Calculated  values  of  the  terms  of  equation  (3.1) 
indicate  that  entrainment  is  the  dominant  term  most  of  the 
time  during  deepening.  In  this  example  the  entrainment  rate 
has  bean  enhanced  by  conditions  set-up  by  the  interior 
motion.  Equation  (3.6)  indicates  that  entrainment  rate.  We, 
is  inversely  proportional  to  mixed  layer  depth,  h,  and 
entrainment  zone  temperature  jump,  <iT.  Compare  Figure  A. 2 
to  Figure  A.1  (the  case  with  no  interior  motion)  anc  A.  3 
(the  270c  phase  case).  It  can  be  seen  that  at  the  start  of 
the  last  quarter  of  the  year  the  90°  phase  case  has  a  much 
lower  mixed  layer  temperature  than  the  other  two  cases  (the 

reason  for  this  will  be  discussed  in  the  next  section)  . 

This  lower  mixed  layer  temperature  means  that  the  entrain¬ 
ment  zone  temperature  jump  will  be  smaller,  allowing  for  a 
larger  entrainment  rate.  At  the  start  of  the  last  quarter 
of  the  year,  the  differences  in  h  are  not  as  significant  as 
those  in  ^T. 

For  example,  the  values  of  il,  h.  We  and  wk  have 

keen  calculated  for  these  three  cases  for  day  280.  The 

mixed  layer  depth,  h,  is  about  St)  m  for  the  case  without 
interior  motion,  65  m  for  the  90e  phase  and  68  m  for  the 
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270°  phase.  These  values  do  not  greatly  differ.  However, 
the  temperature  jump  for  the  90=  phase  case  is  1.3‘C  which 
is  less  than  half  the  2.9'C  value  for  the  case  with  no 
interior  motion  which  ir.  turn  is  less  than  the  4.7'rc  value 
for  the  270'  case.  The  entrainment  rats  for  the  90-'  phase 
case  is  90  cm/day  compared  to  43  cm/day  for  the  no  interior 
motion  case  and  22  cm/day  for  the  270  c  case.  For  both 
phases  considered,  the  interior  motion  is  at  least  an  order 
of  magnitude  smaller  than  the  entrainment  rate. 
Nevertheless,  at  day  280,  the  deepening  in  the  90 J  phase 
case  is  helped  a  little  by  downwelling  and  in  the  270"  case, 
hindered  a  little  by  upwsiling. 

In  the  last  few  days  of  the  year,  the  increase  in 
mixed  layer  depth  begins  to  take  over,  particularly  for  the 
90°  case.  This  inhibits  further,  entrainment  and  the  mixed 
layer  then  rapidly  shallows  to  rhe  cyclical  steady-state 
value  of  75  m  (Figure  A. 9). 

In  the  270"  phase  case  (Figure  A.  3)  there  is  down- 
welling  during  the  middle  half  of  th=  year,  which  increases 
the  deepening  of  the  mixed  layer.  By  the  start  of  the  last 
quarter  cf  the  year,  the  mixed  layer  is  warmer  than  in  the 
case  with  no  interior  motion  or  m  the  90 :  phase  case.  This 
leads  to  a  larger  value  of  ai  and  a  reduced  entrainment 
rate . 

Comparing  Figures  A. 8  and  A.  9  it  can  be  seen  that 
mid-year  upwalling  is  associated  with  a  reduced  rate  of 
deepening  and  mid-year  downwelling  is  associated  with  an 
increased  deepening  rate.  Looking  at  Figure  A. 10,  it  can  be 
seen  that  the  actual  decreased  or  increased  values  of  mixed 
layer  depth  lag  the  start  of  upw<3iling  or  downwelling 
(respectively)  by  90°.  The  time  rate  of  change  of  mixed 
layer  depth  tracks  along  with  the  interior  motion  in 
accordance  with  equation  (3.1).  However,  actual  mixed  layer 
depth  is  the  integral  of  this  and  integrals  of  periodic 
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functions  lag  them  by  90°  .  Also  note  that  the  strong 
enhancement  of  and  strong  hindrance  to  deepening  at  the  end 
of  the  year  (Figure  A.  10)  lags  tae  strong  decrease  and 
strong  increase  in  mixed  layer  temperature  (Figure  A.  12)  by 
about  90". 

2.  Effect  on  “ i xed  Layer  Temperature 

Consider  now  the  mixed  layer  heat  equation  (3.2). 
The  net  surface  heat  flux  is  prescrioed  and  unaffected  by 
the  interior  motion.  The  entrainment  heat  flux,  giver,  by 
equation  (3.6),  is  inversely  proportional  to  mixed  layer 
depth.  The  other  terms  in  equation  (3.6)  depend  primarily 
on  surface  boundary  conditions.  So,  the  first  term  in  equa¬ 
tion  (3.2)  is  inversely  proportional  to  mixed  layer  depth, 
h,  and  the  second,  inversely  proportional  to  the  square  of 
h.  This  implies  that  an  increase  in  h  will  reduce  the 
magnitude  of  mixed  layer  temperature  change;  a  decrease  in  h 
will  increase  the  change  in  TMl_.  This  is  because,  wi-h  a 
deeper  mixed  layer,  a  given  amount  of  heat  added  to  or 
removed  from  tha  mixed  layer  must  be  spread  over  a  larger 
mass  cf  water. 

This  effect  is  illustrated  ir.  the  examples  shown  in 
Figures  A. 2  and  A. 3.  In  Figure  A. 2,  the  mixed  layer  is 
shallower  than  in  Figure  A.  3  from  about  lav  130  to  about  day 
300.  This  resulted  from  the  mid-year  upweiiing  for  the  90J 
phase  case  and  the  aid-year  downwalling  for  the  270-  case. 
The  shallower  mixed  layer  in  Figure  A. 2  allows  for  a  much 
more  rapid  temperature  drop  than  occurs  with  the  deeper 
mixed  layer  of  Figure  A.  3.  This  can  be  seen  in  the  enhanced 
gradients  in  Figure  A. 2.  After  day  300  ,  the  rapid  deepening 
of  the  mixed  layer  for  the  90e  phase  case  and  slow  deepening 
for  the  270°  case  reverse  the  situation. 

To  see  the  effect  of  the  interior  motion  for  all  of 
the  phases,  the  T contours  should  not  be  compared  directly 
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to  the  interior  motion  contours  but  rather  to  the  contours 
of  the  mixed  layer  depth  difference,  Jv  h  (defined  by  equa¬ 
tion  (3.9)).  Between  the  0C  and  133"  phases,  after  day  180 
there  is  a  region  in  Figure  A. 10  where  the  mixed  layer  is 
shallower  with  the  interior  motion  than  without.  This 
causes  the  mixed  layer  temperature  decrease  (the  second  half 
of  the  year  is  the  cooling  season)  to  be  greater,  resulting 
in  the  negative  values  of  ^iTML  in  the  last  quarter  year  at 
these  phases  (Figure  A. 12).  Between  the  ISO01  and  360: 

phases  the  opposite  occurs. 

At  the  end  of  the  year  the  mixed  layer  depth  returns 
to  75  m.  Mixed  layer  temperature  does  not  return  to  its 
initial  value.  The  second  year  T„L  data  is  plotted  in 

Figures  A. 13  and  A. 14.  The  reduced  temperatures  for  phases 
0°  to  180"  and  increased  temperatures  for  phases  180°  to 

360 c  are  continued  into  the  second  year.  The  mixed  layer 
depth  follows  the  same  pattern  ir.  the  second  year  as  it  did 
in  the  first.  This  causes  the  pattern  begun  in  the 

first  year  to  be  enhanced  in  the  second  year.  Continued 
pumping  of  the  ocean  by  this  period  interior  motion  enhances 
this  structure.  If  permitted,  it  would  eventually  cause  the 
upper  200  m  to  become  isothermal  (ir.  the  region  of  the  90° 
phase  at  the  end  of  the  year) ,  followed  by  deepening  beyond 
the  limits  cf  what  the  model  can  handle.  It  would  not  be 
realistic,  however,  to  have  this  single  period  interior 
motion  continue  for  the  several  years  necessary  for  this  to 
happen.  In  reality,  the  interior  motion  would  probably  be 
damped  by  interactions  with  the  boundary  layer  mixing 
processes.  There  should  also  be  more  than  one  period 
present.  Furthermore,  surface  heat  flux  shou:  1  be  influ¬ 
enced  by  sea  surface  temperature.  These  three  possibilities 
open  up  avenues  for  further  research. 
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cnange  ir. 


3.  Effect  on  Heat  Content 


Equation  (3.3)  indicates  that 
content  depends  on  net  surface  heat  flux 
the  ocean  by  the  interior  notion.  The  net 
is  prescribed  and  is  independent  of  the 
Subtracting  equation  (3.3)  for  the  cas 
notion  from  the  same  equation  with  interio 


the  heat 
and  on  pumping  cf 
surface  heat  flux 
interior  motion, 
e  of  r.o  interior 
r  motion  gives 
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(3.  13) 


where  AH  is  defined  by  equation  (3.11).  Equation  (3.13) 
indicates  that  upwelling  will  cause  aH  to  decrease  and 
dowr.welling  will  cause  ah  to  increase.  This  can  be 
observed  by  comparing  Figures  A.  8  and  A.  16.  The  input  to 
heat  content  due  to  net  surface  heat  flux  alone  causes  H  to 
increase  during  the  first  half  of  the  year  and  decrease 
during  the  second  half  as  shewn  in  Pigure  A. 6.  The  combina¬ 
tion  of  surface  heating  with  the  effect  of  the  interior 
motion  produces  the  saddle-shaped  curve  shewn  in  Figure 
A.  15. 

The  effect  of  different  phases  of  the  interior 
motion  can  be  demonstrated  by  the  following  two  simple  exam¬ 
ples  (Pigures  3.2  and  3.3)  .  Ho  mixing  is  permitted  and  a 
two-layer  ocean  is  assumed  for  simplicity. 

In  the  first  example  shown  in  Figure  3.2  there  is 
initial  downwelling  followed  by  heating,  then  upwelling 
followed  by  cooling  and  final  downwelling  to  the  initial 
mixed  layer  depth.  When  downwelling  precedes  heating,  the 
heat  is  distributed  ever  a  thicker  mixed  layer.  A  fixed 
amount  of  heat  will  produce  a  smaller  increase  in  T  than 
would  have  resulted  without  the  downwelling.  If  this  is 
followed  by  upwelling  before  cooling,  the  thinner  mixed 
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layer  will  reach  a  lower  temperature  for  a  given  amount  of 
heat  removed.  At  the  and  of  the  complete  cycle  the  mixed 
layer  temperature  will  be  lower.  Recall  that  the  heat 
content  is  defined  by 

H  -  j  (  T  -  T>  )  a  *  . 

-  D 

Therefore  the  lower  mixed  layer  temperature  produces  a  lower 
heat  content. 

Opweiling  preceding  the  heating  season  distributes 
the  heat  ever  a  thinner  mixed  layer.  This  allows  a  greater 
increase  in  T  than  would  otherwise  occur.  If  this  is 

followed  by  dovnwelling  preceding  the  cooling  season,  a 
fixed  amount  of  heat  removed  will  produce  a  smaller  drop  in 
T.  The  net  effect  will  be  a  higher  temperature  in  the  mixed 
layer  and  a  higher  heat  content. 

Because  of  the  asymmetry,  as  illustrated  in  these 
examples,  heat  is  net  conserved  for  a  particular  phase  of 
the  interior  motion.  This  can  be  seen  in  Figures  A. 15  and 
A. 16.  At  some  phases  the  year-end  heat  content  is  higher 
than  it  was  at  the  beginning  of  the  year,  at  other  phases, 
lower.  However,  if  the  heat  content  is  averaged  over  all  of 
the  phases,  the  year-end  value  eguais  the  initial  value. 
The  average  heat  content  for  sixteen  equally  spaced  phases 
from  0C  to  363°  was  tabulated.  The  phase-averaged  value  of 
H  was  0 . 07 " C  m  less  at  the  end  of  the  year  than  it  was 
initially.  This  small  error  can  be  attributed  to  the  coarse 
resolution  of  averaging  only  sixteen  phases  over  the  entire 
360". 


4.  2.2  Potent ial  Energy 

Time  rate  of  change  of  potential  energy  is  giver,  by 
equation  (3.4) .  From  this  equation  it  can  be  seen  that 
surface  heating  tends  to  increase  potential  energy,  surface 
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cooling  tends  to  decrease  it.  The  net  surface  heat  flux  is 
independent  of  the  interior  motion  so  the  effect  on  :\V  due 
to  this  term  depends  on  ^h.  If  there  is  an  increase  it 
mixed  layer  depth  during  heating,  the  buoyancy  added  will  be 
mixed  deeper,  increasing  potential  energy.  If  there  is  an 
increase  in  mixed  layer  depth  during  cooling  the  effect  will 
be  opposite. 

The  entrainment  heat  flux  term  has  a  positive  sign 
ir.  equation  (3.4).  Increased  mixing  moves  buoyant  water 
down  and  less  buoyant  water  upward,  increasing  potential 
energy. 

Com  par  is  on  of  values  of  the  terms  of  equation  (3.4) 
indicates  that  the  third  term  on  the  right  is  generally 
dominant.  The  exception  is  when  the  interior  motion  goes  to 
zero,  twice  a  year.  The  first  term  of  equation  (3.4)  will 
not  differ  greatly  with  and  without  interior  motion  unless 
there  is  a  large  difference  in  mixed  layer  depth.  Recall 
that  entrainment  heat  flux  is  nearly  inversely  proportional 
to  mixed  layer  depth  (equation  (3.6)).  Thus,  the  second 
term  will  net  vary  significantly.  So,  the  potential  energy 
closely  follows  the  interior  motion.  This  can  be  seer,  by 
comparing  Figures  A. 8  and  A.  18.  During  downweiling, 
potential  energy  increases,  during  upwelling  it  decreases. 

The  physical  reasoning  behind  the  effect  cf  the 
third  term  of  equation  (3.4),  the  pumping  term,  is  that 
downweiling  pushes  buoyant  surface  water  downward,  whereas, 
upwelling  is  in  the  same  direction  as  the  buoyant  forces. 

5.  Summary  qf  Bespon se  to  One -Year  Period  Interior 

get  lot) 

During  the  first  quarter  year  there  is  shallowing 
and  the  interior  motion  has  little  effect  on  either  mixed 
layer  depth  or  temperature.  Bid-year  upwelling  reduces 
deepening  and  mid-year  downweiling  enhances  deepening.  The 
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shallower  mixed  layer  resulting  from  upwelling  will  cool 
faster;  the  deeper  mixed  layer  resulting  from  dowr.wellir.g 
will  cool  slower.  A  cooler  mixed  layer  has  a  reduced  jT 
allowing  for  a  larger  entrainment  rate  and  more  rapid 
deepening.  A  warmer  mixed  layer  deepens  more  slowly. 

Potential  energy  and  heat  are  both  directly  affected 
by  the  interior  metier.  Beth  are  increased  by  dewnweliing 
and  decreased  by  upwelling. 

Before  examining  the  mixing  response  tc  higher 
frequencies,  an  observation  can  be  made  regarding  hcri cental 
structure  produced  by  the  interior  motion.  As  noted  in 

chapter  one,  phase  may  be  related  to  horizontal  displacement 
by  substituting 


Jlf  X 


(3.  14) 


into  equation  {2.20)  to  obtain 
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The  wavelength,  L,  might  be  on  the  order  of  300  km,  large 
enough  for  the  one-dimensio  nal  assumption  to  be  valid. 

Figures  A.  12  and  A.  14  show  the  development  of  hori¬ 
zontal  bands  of  positive  and  negative  mixed  layer  tempera¬ 
ture  anomalies  which  would  result  from  a  one-year  period 
interior  wave.  Because  temperature  is  assumed  to  be  verti¬ 
cally  homogeneous  in  the  mixed  layer,  TMl  also  represents 
sea  surface  temperature.  Figures  A.  11  and  A. 13  show  a  time 
sequence  of  the  horizontal  sea  surface  temperature  field 
resulting  from  the  interior  motion,  the  heat  cycle  and 
mixing  processes.  The  entire  surface  warms  during  the 
heating  season  and  cools  during  the  cooling  season.  During 
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the  first  half  of  the  first  year,  there  is  no  significant 
difference  in  sea  surface  temperature  along  the  horizontal 
axis.  Rfter  the  first  half  year,  the  region  where  there  is 
upwelling  at  the  start  of  the  cooling  season  becomes  and 
remains  cooler  than  the  region  where  rhere  is  dcwnwelling  at 
the  start  cf  the  cooling  season. 

0.  EFFECT  OF  SHOBTEB  PEBIOD  IHTEHIOH  HOTIOM 

Comparison  of  Figures  A. 19  and  A.  21  shows  that  as  before 
downwelling  increases  deepening  and  upwelling  reduces  deep¬ 
ening.  This  car.  also  be  observed  by  comparing  Figures  A. 30 
and  A. 32  and  by  comparing  Figures  A. 39  and  A. 41.  The  mere 
rapid  variation  in  mixed  layer  depth  dees  not  yield  the 
strong  effect  on  temperature  four.!  for  the  one-year  period 
interior  motion.  This  can  be  seen  in  the  weaker  gradients 
of  cf  Figures  A. 23,  A. 34  and  A. 43  compared  to  that  of 
A.  12.  Hence,  the  higher  frequency  interior  motion  does  not 
have  as  strong  an  effect  on  the  mixed  layer  as  the  lower 
frequency  motion  does. 

The  second  year  temperature  data  was  plotted  for  the 
half-year  period  interior  motion  in  Figures  A.  24  and  A. 25. 
Notice  that  the  same  structure  seen  in  the  first  year  reap¬ 
pears  in  the  second,  but  the  cold  spots  have  not  grown 
significantly  colder  nor  the  warm  spots  warmer  as  in  the 
case  of  the  one-year  period  motion. 

As  before,  the  heat  and  potential  energy  differences, 
AH  and  AP,  track  along  with  the  inferior  motion  as  seen  in 
Figures  A. 27,  A. 29,  A. 36,  A. 38,  A. 45  and  A. 47.  The 
resulting  contours  of  H  and  P  appear  more  complex.  This  is 
because  the  values  cf  H  and  P  without  the  interior  morion 
follow  a  cne  year  cycle  {Figures  A. 6  and  A. 7)  and  the  higher 
frequency  cycles  due  to  the  pumping  by  the  interior  motion 
are  superimposed  on  the  one  year  cycle. 
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Fcr  the  cne-year  period  interior  notion,  upwelling  lasts 
long  enough  to  enhance  cooling  significantly  (at  phases 
around  90°)  and  downwelling  lasts  long  enough  to  inhibit 
cooling  significantly  (at  phases  around  2701')  .  The  shorter 
periods  dc  not  maintain  upwelling  and  downwelling  long 
enough  to  produce  an  effect  this  strong. 

E.  ASYHHETBY  IE  ENH1NCED  A  HD  REDUCED  SIZING 

In  Figures  A. 48  and  A. 49  the  values  of  peak  enhanced 
cooling  and  peak  reduced  mixed  layer  cooling  are  plotted. 
Both  first  and  second  year  values  are  considered  (in  Figures 
A. 48  and  A. 49,  respectively).  Negative  values  of  _\T..„, 
indicate  enhanced  coding  due  to  the  interior  motion.  The 
maximum  absolute  value  of  all  of  the  negative  ATMC  data 
points  was  determined  to  indicate  the  peak  enhanced  cooling. 
Positive  values  of  indicate  reduced  cooling.  The 

maximum  value  of  all  of  tha  positive  aTml  data  points  indi¬ 
cates  the  peak  reduced  cooling.  The  positive  values  dc  not 
generally  indicate  enhanced  heating,  because  the  wave  has 
very  little  effect  on  the  mixed  layer  temperature  during  the 
heating  season. 

The  first  obvious  observation  that  can  be  made  from 
Figures  A. 48  and  A. 49  is  that  the  largest  peak  anomalies 
occur  with  the  one-year  period  interior  motion.  For  the 
one-year  period  interior  motion  tha  peak  enhanced  and  peak 
reduced  mixed  layer  cooling  is  2.04  and  2.16°C,  respec¬ 
tively.  For  the  one-half  year  period  motion  •‘these  values 
are  1.35  and  1.15°C,  respectively.  Also,  for  each  of  the 
periods  studied,  the  difference  between  the  peak  enhanced 
cooling  and  the  peak  reduced  cooling  is  greater  in  the 
second  year  than  in  the  first.  All  of  the  periods  studied 
except  the  one-year  period  have  greater  peak  enhanced 
cooling  than  peak  reduced  cooling. 
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The  phase-averaged  values  of  the  variables  Ah,  4T__, 
AH,  and  A?  were  computed  for  sixteen  equally  spaced  phases 
from  0*  to  360"'  .  The  phase-averaged  values  at  the  end  of 
the  first  and  the  end  of  the  second  model  years  are  listed 

in  Table  I.  These  values  indicate  the  net  effect  of  the 

interior  motion  after  each  year  of  interacting  with  the 
mixing  processes.  Consistent  with  the  results  shown  in 
Figures  A. 48  and  A. 49,  the  results  in  Table  I  indicate  that 
the  one-year  period  interior  motion  has  produced  a  net 
reduction  of  mixed  layer  cooling.  The  higher  frequency 
interior  motion  has  produced  a  net  enhancement  of  mixed 
layer  cooling.  For  all  of  the  periods  studied  there  has 
been  a  net  reduction  of  the  mixed  layer  depth.  Heat  is 
conserved,  hence  the  phase-averaged  value  of  AH  is  zero  at 
the  end  of  each  year.  The  effect  of  the  one-year  period 

interior  motion  was  to  decrease  the  potential  energy.  The 

higher  frequency  interior  motion  produced  a  net  increase  in 
potential  energy. 

A  net  increase  (decrease)  in  potential  energy  indicates 
that  the  interior  motion  has  enhanced  (reduced)  the  overall 
mixing.  The  amount  of  deepening  required  to  produce  ar. 
equivalent  increase  in  potential  energy  is  instructive. 
This  can  easily  be  calculated  for  the  simplified  example  of 
a  two-layer  ocean,  illustrated  in  Figure  3.4  with  initial 
mixed-layer  temperature,  T0  and  depth,  h..  This  mixed  layer 
is  deepened  to  h,  with  new  mixed  layer  temperature,  T,  . 
Heat  is  conserved,  so,  using  the  definition  for  the  heat 
variable,  H,  given  in  equation  (2.32)  : 

(r.  -td  )w,«  (t0  -  t B  n,  <3-16> 

Using  equation  (2.31),  the  change  in  the  potential  energy  is 
given  by 
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P  i  -  -  C  1  i  ~  I  n  )  ^  i 


C  Tc  -  T0  '  U  r 


Substituting  aquation  (3.16)  into  aquation  (3.17)  yields 


Po  =  -r(Tc  '  'a)  '  ‘'O 


(3.  18) 


The  term  T„-Tn  is  just  the  initial  entrainment  zone  tempera¬ 
ture  jump.  The  increase  in  potential  energy  due  to  deep¬ 
ening  with  heat  conserved  is  diractly  proportional  to  ths 
increase  in  mixed  layer  depth  for  the  two-layer  example. 


D  l|  _ 


initial  profile 


after  mixing 


Figure  3.4  Two-Layer  Ocean  Example. 


Consider  now  the  phase-averaged  potential  energy  anoma¬ 
lies  at  the  end  of  the  first  year,  shown  in  Table  I. 
Suppose  the  initial  mixed  layer  depth  is  75  m  and  the 
initial  temperature  jump  is  4  °C.  The  net  reduction  of 
potential  energy  for  the  one-year  period  interior  motion 
would  then  be  equivalent  tea  reduced  mixed  layer  depth  of 


5 


0.73  n.  The  net  increase  of  potential  energy  for  the  oro- 
half#  one-third#  and  one-quarter  year  period  interior  action 
is  equivalent  to  an  increase  of  mixed  layer  depth  by  1.6, 
1.8  and  1.8  a#  respectively. 

The  aeount  of  dcwnwelling  and  upwelling  are  the  same 
over  the  whole  year  and  for  all  phases.  Yet,  the  net  reduc¬ 
tion  of  mixing  and  nixed  layer  cooling  does  not  equal  the 
net  enhancement  of  the  same.  The  net  affect  of  the  one-year 
period  interior  motion  is  to  hinder  mixing.  The  net  effect 
of  the  shorter  period  interior  motion  :s  to  enhance  mixing. 


TABLE  I 

Phase-averaged  Bean  Anomalies 
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Da  y 

&h 

AH  | 

A? 
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1 
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0, 
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0. 
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0. 
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IV.  COBCIOSIONS  iNS  R ECO H HEIDATI01S 


A.  COHCIOSIOBS 

Seasonal-scale  interior  motion  has  a  significant  effect 
on  boundary  layer  mixing  processes.  This  effect  is  most 
pronounced  when  the  period  of  the  interior  motion  is  one 
year.  For  this  period,  upwelling  between  the  heating  and 
cooling  seasons  can  strongly  enhance  mixed  layer  cooling  and 
subsequent  deepening.  Downwelling  can  strongly  reduce  the 
same.  This  effect  may  produce  a  significant,  observable 
horizontal  pattern  in  sea  surface  temperature.  Although 
there  is  negligible  effsct  on  sea  surface  temperature  and 
mixed  layer  depth  during  the  beginning  of  the  heating 
season,  the  interior  motion  dees  affect  tne  potential  energy 
ar.d  heat  content  at  this  time.  Throughout  the  year  the 
potential  energy  and  heat  are  directly  related  to  the 
pumping  by  the  interior  motion.  Upwelling  decreases  and 
downwelling  increases  both  potential  energy  and  heat 
content. 

Higher  frequency  interior  motion  does  not  have  as  strong 
an  effsct.  The  cycling  between  upwelling  and  downwelling 
during  the  cooling  season  prevents  the  strongly  enhanced 
mixed  layer  cooling  from  occurring. 

The  effect  of  the  interior  motion  is  asymmetric. 
Generally,  the  enhancement  of  mixing  and  mixed  layer  cooling 
exceeds  the  reduction  of  the  same.  The  exception  to  this 
occurs  for  the  one-yaar  period  intarior  motion. 

B.  HECOHHE1DATIOMS 

One  aspect  of  this  research  which  departs  from  reality 
is  in  prescribing  a  single  wave  for  the  interior  motion. 
This  could  be  changed  quit9  easily  by  prescribing 
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w. 


-  *»)■ 


(4.1) 


The  amplitudes,  and  periods,  T* ,  could  be  based  on  the 

first  few  dominant  frequency  components  observed  in  actual 
wave  spectra. 

Another  way  that  reality  could  be  more  closely 
approached  is  by  allowing  the  net  surface  heat  flux  to  vary 
depending  cn  the  sea  surface  temperature.  This  should 
dampen  the  strongly  enhanced  or  reduced  cooling  of  the  mixed 
layer  temperature  seen  with  the  one-year  period  motion. 

An  area  open  to  further  research  is  to  use  a  finite 
depth  ocean  and  an  appropriate  vertical  wave  structure.  An 
appropriate  wave  structure  for  planetary  waves  is  giver,  by 
Hillmott  and  Mysak  (1980)  .  In  their  work  the  vertical 
motion  is  linear  in  the  mixed  layer: 


i)1  w 

<)  i  “ 


C  , 


(*.2) 


but,  below  the  mixed  layer  the  vertical  motion  is  given  by 


N  i- 
O 


Vv 


o 


(9-3) 


where  N2  is  the  3runt-Vaisaia  frequency  and  c„  represents 
the  speed  of  propagating  long  waves  in  the  non-rctatir.g  case 
for  mode  number  n.  The  solutions  of  equations  (4.2)  and 
(4.3)  must  be  matched  at  the  base  of  the  mixed  layer,  z=-h, 
to  ensure  continuity  cf  the  vertical  velocity.  The  boundary 
conditions  are  rigid  lid  and  no  normal  flow  through  the 
ocean  bottom,  so,  w*0  at  z=0  and  at  the  ocean  bottom.  At 
each  time  step  the  mixed  layer  depth  is  recalculated  and  so 
the  vertical  wave  structure  problem  would  also  have  to  be 
solved  for  the  new  w  (z)  .  The  combination  of  this  sort  of 
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vertical  wav®  structure  with  the  Garwood  model  for  boundary 
layer  mixing  processes  could  provide  an  interesting  area  for 
future  study. 

In  this  study  the  interior  motion  was  prescribed  and 
acted  on  the  mixed  layer.  The  interior  motion  was  not 
allowed  to  be  changed  by  the  mixed  layer  processes.  If,  for 
example,  the  wave  motion  enhances  mixing,  the  increase  in 
potential  energy  should  be  balanced  by  a  decrease  in  wave 
energy.  Although  the  frequencies  considered  here  are  much 
lower  than  those  in  the  work  by  Bell  (1978),  the  processes 
involved  could  bs  sixilar. 


0  30  60  90  120  150  180  210  240  270  300  33C  360 

TIME  DRY5 

Figure  A. 2  Temperature ,  One-Tsar  Period,  90°  Phase. 
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Figure  A. 3  Temperature,  one-lear  Period,  270°  Phase 


56 


o 


o 


Figure  1.5  sized  Layer  Tenperatura,  lo  Interior  Botion 
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Figure  1.6  Heat,  Ho  Interior  notion. 
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T)MC  P*X5 


Figure  1.8  One- Tear  Period  Interior  Hotion. 
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Figure  A. 13  TMl_  (C) ,  One-Tear  Period,  Second  Tear 
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Figure  1.14  One-Tear  Period,  second  Tear. 
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T  WE  DRYS 


Figure  A.  15  H  (100‘C  ■) ,  One-rear  Period  Interior  lotion. 
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Figure  A.  16  AH  (100cC  a),  one-fear  Period  interior  lotion. 
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Figure  *.18  AP  (10*X  one-lea r  Period  interior  notion. 


d 


Figure  A. 19  1/2-Year  Period  Interior  Hotion. 
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f*gar«  1.27  &H  (100'c  ■)  ,  1/2-I«ac  Period  interior  action. 
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Figure  k.  28  P  (10*  "C  ■*)  ,  1/2- Tear  Period  Interior  notion. 
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Pigure  i.29  ^P  (10* °C  ■*),  1/2-Iear  Period  interior  notion, 


Figure  A. 32  Ah  (meters) ,  1/3-Yee.r  Period  Interior  Hot  ion. 
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TIME  DAYS 

Figure  a. 35  H  (100  *C  a),  1/3-Year  Period  Interior  notion 
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figure  a. 44  H  (1Q0“C  ■) ,  1/4-Tear  Period  Iaterior  notion. 


Figure  A. 45  & H  (100°c  ■)  ,  1/4-Tear  Period  Interior  notion 


Figure  A. 48  Peak  Enhanced  and  Reduced  Cooling  Values,  Tear  1 


Figure  A. 49  Peak  Enhanced  and  Reduced  Cooling  Values,  Tear  2 
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